Few particles confined in an open pore 
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Abstract 

This paper presents a modified grand canonical ensemble which provides a new simple and efficient 
scheme to study few-body fluid- like inhomogeneous systems under confinement. The new formalism 
is implemented to investigate the exact thermodynamic properties of a hard sphere (HS) fluid-like 

r^ \ system with up to three particles confined in a spherical cavity. In addition, the partition function 

O ' 

CJ ' of this system was used to analyze the surface thermodynamic properties of the many-HS system 
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and derive both, the exact curvature dependence of the surface tension and adsorption, in powers 
of the density. We use these results to derive for this system the dependence of the fluid-substrate 
Tolman length up to first order in density. 



INTRODUCTION 

The thermodynamic and statistical mechanical properties of a fluid composed of many 
particles confined in small pores have been extensively investigated from both a theoretical 
and experimental point of view l|-l3|. In these systems the fluid fllls the cavity inhomoge- 
neously and under extreme conditions, the conflnement is responsible for the dimensional 
crossover of the particles (^. On the contrary, the knowledge about the statistical mechan- 
ical properties of fluid-like systems composed of few bodies is rather poor. Notoriously, 
few-body systems have played a distinguished role in the development of various flelds of 
physical science including classical, relativistic and quantum mechanics. The fundamental 
advantage of the few-body systems over the many-body ones obviously resides in that only 
the former are (in most cases) analytically tractable. In this sense we claim that, the study 
of few-body systems in the framework of the statistical mechanics of inhomogeneous fluids 
is interesting because it provides exact solutions in a subject where there are extremely few 
of such results, thus helping to better understand the properties of many-body systems. 

The grand canonical ensemble (GCE), in which the volume, temperature and chemical 
potential are flxed, is certainly the most used scheme in the formulation of theories of fluids. 
In particular, the density functional theory (DFT) formulated in GCE, is one of the most 
successful approaches to properly describe the properties of inhomogeneous fluids. Indeed, 
DFT calculations make it possible to analyze the adsorption of fluids on substrates and in 
porous matrices as well as study the interfacial phenomena. Additionally, other theories as 
the scaled particle theory (SPT) are formulated in the GCE. In this context the hard spheres 
(HS) systems play an important role in the development of DFT and SPT approaches and 
in the perturbation theories of fluids |5|. Moreover, HS systems are of particular interest as 
they constitute a simplifled model for both simple fluids and colloidal particles (6|. 

In this work we present a novel formulation that makes it possible the study of few-body 
inhomogeneous fluids using a GCE scheme in which the maximum number of particles is 
considered as a parameter of the system. We employ the proposed GCE to perform an 
exact study of the HS system with at most three particles conflned in a spherical cavity. In 
addition, on the basis of the new results, we analyze the surface thermodynamic properties 
of the many-HS system spherically conflned. We present our result on the exact curvature 
dependence of surface (or boundary) tension and adsorption, obtained in powers of the 



density. The formulation proposed here is found to be a valuable tool for describing highly 
confined fluids constituted by few bodies. Moreover, it will potentially serve to improve the 
DFT and SPT descriptions associated with the surface-related thermodynamic properties 
and their curvature dependence. 

PARTITION FUNCTION OF THE FEW BODY OPEN SYSTEM 

Consider a one-component fluid containing a non-fixed number of particles which evolve 
inside a region of the space that we will refer to as B. Such a region is determined by an 
external potential (f{r) that takes finite values in B and diverges outside it. The boundary 
of B is assumed to be at a constant and uniform temperature T. In these circumstances 
we expect the fluid will attain an equilibrium state in which its physical properties either 
remain constant or fluctuate around a fixed value. Furthermore, for an inhomogeneous fluid 
system in an equilibrium state we assume the following standard hypothesis: i) if the exact 
expression for the grand canonical partition function (GCPF) S is known, then the exact 
grand potential function is given by 

Q = -/3~HnE, (1) 

where /3 = {kBT)~ is the inverse temperature (with ks and T being the Boltzmann's 

constant and the temperature, respectively); and ii) a given thermodynamic property of the 

system, namely X, is identified with its mean (Gibbsian) ensemble value {X). Additionally, 

it is expected that {X) to be equal to its time averaged value over an interval r, Xr = 

- J X{t)dt (at least for a value of r large enough). Typically, X may be the energy E, the 

pressure p, the number of particles A^, etc. Standard statistical mechanical demonstrations 

along with usual thermodynamic relations are necessary to conipute from the grand potential 

and its derivatives all of the thermodynamic quantities X p, l7|. 

In order to study few-body inhomogeneous systems, we now consider a container B with 

at most M. particles, which in what follows will be a fixed parameter. The GCPF of this 

system is given by 

M 

S^ = ^I,^^Z,, (2) 

i=o 

where M. is the cutoff value for the maximum number of particles. The factor Ij is the 

so-called indistinguishably factor, which is either equal to I^ = 1/j! for indistinguishable 



particles or Ij = 1 for distinguishable particles. In addition, z = A~^ exp /3/i is the activity 
being A = h/(2TTmkBTY^'^ the thermal de Broglie wavelength, while /i, m and h are the 
chemical potential, the mass of each particle and the Planck's constant, respectively. Finally, 
The configuration integral (CI) of a system with j particles is 

Z, = J ...fl[eil[e,n,d^r, (3) 

being Zq = 1. Here, ei = exp [— /3(y9(r/)] and eim = exp [—f3(f>{rim)] are the Boltzmann's 
factors corresponding to the external and pair interactions, respectively. Note that, the 
integration domain in Eq. ([3]) is the complete space due to the spatial confinement of the 
particles in the region B is considered in the e/ terms. 

In view of the fact that S_a4 is the GCPF of our system, we can replace in Eq. ([1]) S 
with H_A4 to define Q. The exact knowledge of Zj for 1 < j < Ai along with Eq. ([1]) and 
Eq. ([2]) provide the exact statistical-thermodynamic properties of the system. In Eq. (^, 
the expression for S_a4 is an incomplete (or restricted) version of the usual GCPF obtained 
from S = lim S_a4. A different restricted GCE was introduced by Yang and Lee [8J in their 
study of condensation. It was used later by Woods et al. (9| to study the adsorption of 
fluids in cavities. We emphasize that in their analysis Ai represents the maximum number 
of particles that fit in B, instead of an externally imposed parameter as we have assumed in 
our formalism. An interesting property of S_a4, Eq. ([2]), is that for a cavity B with a fixed 
size and temperature, lim Sx oc IjZ^Zj with j the maximum number of particles that holds 

z—>oo 

in B (under the constraint 1 < j < Ai). 

THE FEW HS SYSTEM IN A SPHERICAL CONFINEMENT 

Few-body fiuid-like systems of HS with one or two particles confined in cavities with 
different geometries, including the spherical, cuboidal and cylindrical cases, were analytically 



solved in the canonical ensemble in Refs. 



10 



-|l2|. That work made it possible to explore 



on exact grounds the corresponding open systems with A^ = 2. Further^ the three-body HS 
system was exactly solved for the case of a spherical confinement 13|, llJj. In the current 
work we focus on the study of the HS open system with A^ = 3 in a spherical cavity. 
Inhomogeneous systems with a spherical geometry are ubiquitously present in nature in 
the form of bubbles and drops. One important characteristic of this symmetry is that it 



does not place any preferential direction on the space. In spite of the relevance of the 
spherical geometry much efforts have been directed in recent years towards answering still 
open questions about the curvature dependence of the thermodynamic properties in spherical 
systems. 

In our description we consider an external hard-wall potential, which is null in a spherical 
region B and diverges outside it. The Boltzmann's factor is ei = Q{R — ri), where 6(r) is 
the Heaviside function, r; is the distance from the center of the pore to the particle / and R 
is the effective radius of the pore 13|. The volume oi B is V = Zi = ^R^ and A = AuR^ is 



;he surface area of its boundary. By using the expressions for Z2 and Z3 presented in Refs. 
lOl . I13L Il4| and assuming indistinguishable particles, i.e., Ij = 1/j!, we obtain the exact 



expression for S3. 

Regarding the basic geometrical properties of the spherical confinement we find that, after 
fixing the diameter of the particles to be a = 1 to lighten our notation, a cavity with -R < ^ 
can contain at most one particle, while for the cases i? < 4= = 0.57735, R < -1 /| = 0.612372 
and R < -y= = 0.707107, a maximum of two, three and four particles can be fitted in the 
cavity, respectively. These values imply that by evaluating S3 we have not only solved the 
restricted system with A^ = 3 and any value of R, but also we have obtained the exact 
solution of the unrestricted system for /? < -1 /| because in this case S(i?) = '^^(R). 

The thermodynamic fundamental relations for the grand potential are 

n = u -TS-12N , (4) 

dn = -SdT-PwAdR-Ndfi, (5) 

here S is the entropy, dV is equal to AdR, Pw represents the (mean) pressure on the 
wall also called the work-pressure, and dW = PwdV is the total reversible work performed 
by the system on its environment. Given that the system is athermal, it is preferable to 
use temperature-independent quantities like PVL to conduct the analysis of its properties. 
By following this argument, some of the usual thermodynamic and statistical mechanical 
relations for systems under spherical confinement can be written as 

,1 dpVt 



(3Pw = -A' 



dR 

d(5VL 



dz 



(6) 
(7) 
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Table I. Fugacity and bulk-density values used for the curves plotted in Fig. [TJ 
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dz 



(8) 



Here, a^ is the standard deviation in the number of particles which quantifies the spon- 
taneous fluctuation of A^. Since the energy of the system is equal to that of the classical 
ideal gas, i.e., /3f/ = |A^, using Eqs. (jl]) and ([7]) we can also calculate S. At this stage, 
the exact properties obtained for the system are functions of z and R. One important point 
is the fact that the few-body HS open system is in chemical equilibrium with a bulk HS 
fluid. Unfortunately, the exact properties of the bulk HS fluid are not analytically known, 
and thus, to present our results in terms of the bulk density pb (instead of z) we adopt the 
simple and accurate Carnahan-Starling's equation of state, which for z gives |4l. Il5| 

"% - 9r}l + 2,ril 



^cs = Pb* exp 



(9) 



(1 - v^r 

where 77^ = ^PbO''^ is the corresponding packing fraction. 

Fig. [Udisplays the pressure on the wall as a function of the radius of the cavity, considering 
different values of pb. Table |T] summarizes the values of both z and Pb{z) used in plotting 
each curve. These values were derived from Eq. ([9]). The dashed vertical line in Fig. [1] 
corresponds to R = 0.612, that is, the maximum value for which our results are exact even 
for the unrestricted system. From this graphic it is evident that the pressure on the wall has 
a non-monotonous behavior characterized by regions of negative and positive slopes (van der 
Waals loops), related to the mechanic instability of the system. As a reference, we mention 
that in the bulk HS system the fluid (disordered) and solid (ordered) stable phases coexist 
at B P ~ 11.576 with fluid and solid bulk-densities pf,f — 0.943 and pbs — 1.041, respectively 
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17|. In Fig. [2] we plot the relative adsorption in the pore per unit volume as a function 
of Pb, that is Ap = p — Pb where p = N/V. There, the arrow indicates increasing values of 
R. The negative adsorption in the case of the curves corresponding to R =1, 1.2 and 1.5 
is a consequence of the cutoff in the maximum number of particles contained in the pore. 
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Figure 1. (color-online) Pressure on the wall against the radius of the pore obtained for the values 
of the bulk density, pb, listed in Table [I] (p^ increases from bottom to top). The dashed line indicates 
the radius for which a maximum of three particles can be fitted in the cavity. 




Figure 2. (color-online) Relative adsorption in the pore as a function of the bulk density plotted for 
different radii of the pore R =0.2,0.4,0.5, 0.525,0.55,0.575, 0.6,0.61237, 0.65,0.7,0.8, 1,1.2, and 
1.5. The arrow marks the direction of increasing values of R. The asterisk indicates the curve for 
R = 0.61237. 



The same conclusion applies to the negative slope region for R =0.65,0.7 and 0.8. Fig. [3] 
shows the filling of the cavity as a function of p^ for different radii of the pore. In this plot 
the saturation value of A^ corresponding to each value of R coincides with the maximum 
number of particles (at most three) that holds in the pore. As it can be observed, N{pf), R) 
is a non-decreasing function. It is noticeable that for 0.5 < R < 0.61237 the curves change 
their concavity several times with the increase of pi,. This property determines the behavior 
of cffipb) through Eq. ([8]). Fig. H] illustrates the relation between the mean number of 
particles in the cavity and its fluctuation a^^. We find that the curves corresponding to 




Figure 3. (color-online) Mean number of particles vs. bulk density for different radii of the pore. 
The arrow, the asterisk and the values of R were described in Fig. [2j 




Figure 4. (color-online). Standard deviation in A^ as a function of N for different radii of the pore. 
The arrow, the asterisk and the values of R were described in Fig. [2j 



R < 0.5 collapse onto the semi-circle centered at A^ = 0.5. For 0.5 < R < 0.577 and 
1.5 < N < 2 the curves collapse onto the second semi-circle centered at A^ = 1.5, while 
for R > 0.577 and 2.5 < A^ < 3 the curves collapse onto the third semicircle. In addition, 
we have analyzed the case M. = 2 and obtained a picture (the graphic is not shown here) 
similar to that displayed for the case A^ = 3 but with two semicircles centered at A^ = 0.5 
and A^ = 1.5. One important consequence of the graphic showed in Fig. HJis that it can be 
straightforwardly extended to obtain the representation for the case with a value of A^ > 3. 
Basically, it will involve a sequence of Ai semicircles of unit diameter centered at half-integer 
values of A^ along the abscissa axis. In particular, for a given R the curve will start at the 
origin (A^, a^) = (0, 0), and will collapse either onto the semicircle that ends at the value A^ 
equal to the maximum number of particles that fit in this cavity or onto the A^-th semicircle 
in case R is large enough. It is to be noted that, the outmost curve which connects A^ = 
to A^ = A^, can be obtained analytically. 



THE MANY HS SYSTEM IN CONTACT WITH A SPHERICAL WALL 

In what follows we draw attention to the implications that the expression for Q, obtained 
for the system with A^ = 3, has on the surface-thermodynamic properties of the many-body 
system. In particular, we will address on the HS system in contact with a hard spherical 
wall in the low density limit. For each magnitude of the system, our main goal is to find its 
exact power series in ph up to the three bodies term, where the dependence on the curvature 
is included in the coefficients of the series. Our procedure is based on the ideas developed 



in Ref. 



12|, which enable us to express the Zi{R) terms analytically known, as a function 



of the set of measures {V, A, R}. We begin by placing the dividing surface at R and we use 



the same decomposition rule adopted in Ref. |l3| (Eqs. (28) and (29) therein). For the 
unknown Zi terms we assume a generic dependence Zi{V,A,R). As a result, we obtain an 
expression for fi as a function of V, A, R and z. We then fix A^ = 4 in Eq. ([2]) and calculate 
the magnitude X. We expand the expression for X as a power series in the variable ph by 
following standard procedures based on the inversion and composition of the power series. 
Then, we truncate the series to the first term in which Z4 appears. In addition, to obtain the 
planar-surface property we consider the limit of the series as i? — )■ 00. Thus, we determine 
in the planar limit the expressions for the fluid-substrate surface-tension 7 = dQ/dA\v,R 
and the adsorption per unit area T = V {p — pb) /A, 

/37oo = -^pg(l + ^P6)+0(p^), (10) 

roo = ^P^(i-^P.)+o(p^). (11) 

Eq. (fTOj) exactly reproduces the first two terms (the only known up to present) of the 
series expansion of /37oo(p6)- The above Eqs. are consistent with the Gibbs adsorption 
isotherm Foo = —d'joo/dp. Moreover, Eq. (fTOj) and Eq. (ITT]) are in close agreement both 



density p^ ~ 0.4 (see 



18[) and theoretical results 



with simulation results up to a moderate 

obtained using SPT (see Figs. 2 and 3 in 19|). As is well known, the important problem in 



statistical mechanics regarding the dependence of the surface thermodynamic properties on 
the curvature remains unresolved, even for the HS fluid in contact with hard spherical walls. 
Following the procedure described above and writing the magnitude X/X^o as a double 



power series in the variables R ^ and pb, we obtain 

7(i?)/7oo = 1 + <5oo 2R-^ + S^R-'' + 0{R~^) , (12) 

T{R)/T^ = 1 + ^j 2i?-' + ^ki?"' + 0{R-^) . (13) 

We included in the appendix the complete expression for 7(-R) and T{R) in the limit of 
both small p^ and large R. In the latter equations, 2R~^ and R~^ are the mean and Gaus- 
sian curvatures of the spherical cavity, respectively. The coefficient (5oo denotes the radius- 
independent wall- fluid Tolman length. In fact, neither of the coefficients in the latter two 
equations depends on R. We verified that Eq. flT2|) and Eq. flT3l) hold for fluids in contact 
with a hard spherical wall, irrespective of whether the fluid is inside the cavity (i? > 0) or 
outside of a flxed hard sphere [R < 0, also known as an empty cavity in a bulk fluid) |13| . 
To flrst order in density the coefficients 5oo and 5k can be computed from 

5oo = 0.34299 Pb A = -0.05555 + 0. 19933 pb , 
^j = 0.51448 Pb A = -0.05555 + 0.29899 pb . (14) 

No numerical simulation results have been published up to present for these magnitudes. The 
expressions for both (5oo and (5k determined by the above procedure are consistent with those 



recently reported in 13|, derived from an exact calculation for a few-body system treated 



19| evaluated 



19|). From the comparison 



in the framework of the canonical ensemble. In addition, Siderius et al 

5oQ by using six different versions of SPT (see Fig. 5 in Ref. 

between their and our results we find that only the expression for 5oc obtained by using the 

CS-SPT version is compatible with the low density behavior presented here in Eq. (jT 



Indeed, in Ref. jl9| it is shown that the E-SPT version gives a value Soo{pb = 0) 7^ |20| 
while the others SPT versions yield slopes that differ by a factor of about three from the 
exact value derived from Eq. (IT4|) . At this point, we must note that the dependencies of 7 
and r on the curvature [Eos 



which they were obtained 



12 



fT2lfT^ ] strongly depend on the details of the procedure by 
13| as well as on the choice for the location of the dividing 
surface. Further details regarding the present treatment will be discussed elsewhere. It is 
important to emphasize that the thermodynamic limit was not invoked in our calculations 
of the properties of the many-body system. Moreover, although the zqs equation of state 
given in Eq. Q is used in the actual calculations, it does not introduce any approximation 

10 



in Eqs. (llOtiMp due to zqs is exact to the order of the three-body term involved in our 
results. 

In summary, we presented a restricted GCE form, which is specially useful for analyzing 
the statistical mechanical properties of few-body open systems. We used the partition 
function of this restricted ensemble to evaluate, in an exact analytical framework, several 
properties of the HS fluid-like system composed of at most three particles. Additionally, 
we extended our analysis to derive the surface thermodynamic properties of the many-body 
system of HS in contact with a hard spherical wall. By adopting a power series representation 
we reported new expressions for both the surface tension and adsorption, as a function of 
the density and curvature. We use these results to determine the Tolnian length of the 
fluid-substrate system. 

This work was supported by Argentina Grants CONICET-PIP 0546, UBACyT 20020100200156 
and ANPCyT-PICT 2008-0795. 

APPENDIX 



The complete expressions for /37 and F 
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